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Abstract 

We consider the dynamics of a periodic chain of N coupled overdamped particles under 
the influence of noise, in the limit of large N. Each particle is subjected to a bistable 
local potential, to a linear coupling with its nearest neighbours, and to an independent 
source of white noise. For strong coupling (of the order iV^), the system synchronises, 
in the sense that all particles assume almost the same position in their respective local 
potential most of the time. In a previous work, we showed that the transition from 
strong to weak coupling involves a sequence of symmetry-breaking bifurcations of the 
system's stationary configurations. We analysed, for arbitrary N, the behaviour for 
coupling intensities slightly below the synchronisation threshold. Here we describe 
the behaviour for any positive coupling intensity 7 of order N'^, provided the particle 
number N is sufficiently large (as a function of j/N'^). In particular, we determine 
the transition time between synchronised states, as well as the shape of the "critical 
droplet" to leading order in Our techniques involve the control of the exact 

number of periodic orbits of a near-integrable twist map, allowing us to give a detailed 
description of the system's potential landscape, in which the metastable behaviour is 
encoded. 
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1 Introduction 



In this paper, we continue our analysis of the metastable dynamics of a periodic chain of 
coupled bistable elements, initiated in [BFG06aj . In contrast with similar models involv- 
ing discrete on-site variables, or "spins", whose metastable behaviour has been studied 
extensively (see for instance |dH04[ [OV05] ). our model involves continuous local variables, 
and is therefore described by a set of interacting stochastic differential equations. 

The analysis of the metastable dynamics of such a system requires an understanding of 
its A^-dimensional "potential landscape" , in particular the number and location of its local 
minima and saddles of index 1. In |BFG06a] . we showed that the number of stationary 
configurations increases from 3 to 3'^ as the coupling intensity 7 decreases from a critical 
value 7i of order A^^ to 0. This transition from strong to weak coupling involves a sequence 
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of successive symmetry-breaking bifurcations, and we analysed in detail the first of these 
bifurcations, which corresponds to desynchronisation. 

In the present work, we consider in more detail the behaviour for large particle number 
N. In the limit N ^ oo, the system tends to a Ginzburg-Landau stochastic partial dif- 
ferential equation (SPDE), studied for instance in [EH01|, [Rou02j . The Ginzburg-Landau 
SPDE describes in particular the behaviour near bifurcation points of more complicated 
equations, such as the stochastic Swift -Hohenberg equation |BHP05j . For large but finite 
N, it turns out that a technique known as "spatial map" analysis allows us to obtain a 
precise control of the set of stationary points, for values of the coupling well below the 
synchronisation threshold. More precisely, given a strictly positive coupling intensity 7 
of order N'^, there is an integer No{'y/N'^) such that for ah N ^ iVo(7/^^^), 

we know 

precisely the number, location and type of the potential's stationary points. This allows 
us to characterise the transition times and paths between metastable states for all these 
values of 7 and N. 

This paper is organised as follows. Section [2] contains the precise definition of our 
model, and the statement of all results. After introducing the model in Section 12.11 and 
describing general properties of the potential landscape in Section 12.21 we explain the 
heuristics for the limit ^ 00 in Section [2.31 In Section [2.41 we state the detailed results 
on number and location of stationary points for large but finite N, and in Section 12.51 
we present their consequences for the stochastic dynamics. Section [3] contains the proofs 
of these results. The proofs rely on a detailed analysis of the orbits of period iV of a 
near-integrable twist map, which are in one-to-one correspondence with stationary points 
of the potential. Appendix lAl recalls some properties of Jacobi's elliptic functions needed 
in the analysis, while Appendix [B] contains some more technical proofs of results stated in 
Section 13. 61 
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2 Model and Results 
2.1 Definition of the Model 

Our model of interacting bistable systems perturbed by noise is defined by the following 
ingredients: 

• The periodic one-dimensional lattice is given by A = Z /NI^ , where A ^ 2 is the 
number of particles. 

• To each site i S A, we attach a real variable G M , describing the position of the ith 
particle. The configuration space is thus X = M.^. 
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• Each particle feels a local bistable potential, given by 

UiO = lc^-lf, eeM. (2.1) 

The local dynamics thus tends to push the particle towards one of the two stable 
positions ^ = lor^ = — 1. 

• Neighbouring particles in A are coupled via a discretised-Laplacian interaction, of 
intensity 7/2. 

• Each site is coupled to an independent source of noise, of intensity a. The sources of 
noise are described by independent Brownian motions {Bi{t)}t^Q. 

The system is thus described by the following set of coupled stochastic differential 
equations, defining a diffusion on X: 

dx^it) = f{x1 [t]) dt + I - 2x^{t) + xUit)] dt + adBiit) , i G A , (2.2) 

where the local nonlinear drift is given by 

/(e) = -vc/(e) = e-f ■ (2.3) 

For 0" = 0, the system ()2.2p is a gradient system of the form x = — VKy(x), with potential 

V,{x) = U{x,) + ^ 5](x,+i - x,f . (2.4) 

2.2 Potential Landscape and Metastability 

The dynamics of the stochastic system depends essentially on the "potential landscape" 
Vy. As in [BFG06a] . we use the notations 

S = 5(7) ={x£X: VV^{x) = 0} (2.5) 

for the set of stationary points, and 5*^(7) for the set of /c-saddles, that is, stationary points 
with k unstable directions and N — k stable directions. 

Understanding the dynamics for small noise essentially requires knowing the graph 
G = {So,£), in which two vertices x*,y* £ Sq are connected by an edge e G f if and only 
if there is a 1-saddle s £ Si whose unstable manifolds converge to x* and y*. The system 
behaves essentially like a Markovian jump process on Q. The mean transition time from 
X* to y* is of order q^^/'^ , where H is the potential difference between x* and the lowest 
saddle leading to y* (see |FW98] ). 

It is easy to see that S always contains at least the three points 

O = (0,...,0), /± = ±(1,...,1) . (2.6) 

Depending on the value of 7, the origin O can be an A^-saddle, or a A;-saddle for any odd 
k. The points always belong to 5o, in fact we have 

V^{x)>V^{I+)=V^{r) = -^ Vxe A" \ {/-,/+} (2.7) 

for all 7 > 0, so that and represent the most stable configurations of the system. 
The three points O, and /~ are the only stationary points belonging to the diagonal 

V = {x £ X : xi = X2 = ■ ■ ■ = xn} ■ (2.8) 
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N 


X 


Type of symmetry 


AL 


A 


(Xi, . . . ,XL,XL, . . . ,Xi 


, -Xl,.. 


■ , -XL, -XL, -Xl) 




B 


(Xi, . . . ,XL, . . . ,Xi,0, - 


-Xl, . . . 


-XL, -Xi,0) 


4L + 2 


A 


(Xi, . . . , X2,+l, . . . , Xl, - 


-Xl, . . . 


, -Xi+i,...,-Xi) 




B 


(Xi, . . . ,XL,XL . . . ,Xi, 


0, -Xl, 


. . , -XL, -XL, -Xi,0) 


2L + 1 


A 


, . . . , ^L; • • • ; ~ 


-xi,0) 






B 


(Xi, . . . ,XL,XL, . . . ,Xi 


,xq) 





Table 1 . Symmetries of the stationary points bifurcating from the origin at 7 = 71 . The 
situation depends on whether N is odd (in which case we write TV = 2L + 1) or even (in 
which case we write N ~ 4L 01 N — 4:L + 2, depending on the value of N (mod 4)). Points 
labelled A are 1-saddles near the desynchronisation bifurcation at 7 = 71 , those labelled 
B are 2-saddles (for odd N, this is actually a conjecture). More saddles of the same index 
are obtained by applying elements of the symmetry group Gn to A and B. 



On the other hand, being a polynomial of degree 4 in variables, the potential V-y can 
have up to 3'^ stationary points. 

The potential Ky(x), as well as the sets S{'~f) and Skij), are invariant under the trans- 
formation group G = Gn oi order 4:N (4 if N = 2), generated by the following three 
symmetries: 

• the rotation around the diagonal given by i?(xi, . . . ,xn) = {x2, ■ ■ ■ ,xn,xi); 

• the mirror symmetry 5(xi, . . . , x^v) = {xn, ■ ■ ■ , xi); 

• the point symmetry C(xi, . . . , xat) = — (xi, . . . , x^v). 

In |BFG06a] . we proved the following results: 

• There is a critical coupling intensity 

1 — cos(27r/A') 

such that for all 7 ^ 71, the set of stationary points S consists of the three points O 
and only. The graph Q has two vertices connected by a single edge. 

• As 7 decreases below 71, an even number of new stationary points bifurcate from 
the origin. Half of them are 1-saddles, while the others are 2-saddles. These points 
satisfy symmetries as shown in Tabled) The potential difference between and the 
1-saddles behaves like A''(l/4 — (71 — 7)^/6) as 7 /" 71. 

• New bifurcations of the origin occur for 7 = = (1 ~ cos(27rM/A^))^"'^, with 2 ^ 
M ^ N/2, in which saddles of order higher than 2 are created. 

The number of stationary points emerging from the origin at the desynchronisation 
bifurcation at 7 = 71 depends on the parity of N. If N is even, there are exactly 2N new 
points (N saddles of index 1, and N saddles of index 2). If is odd, we were only able to 
prove that the number of new stationary points is a multiple of AN, but formulated the 
conjecture that there are exactly 4A^ stationary points (2A^ saddles of index 1, and 2N 
saddles of index 2). We checked this conjecture numerically for all N up to 101. As we 
shall see in Section 12.41 the conjecture is also true for N sufficiently large. 
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2.3 Heuristics for the Large-A^ Limit 

We want to determine the structure of the set S of stationary points for large particle 
number N, and large coupling intensity 7. For this purpose, we introduce the rescaled 
coupling intensity 

7 27r2 



^71 iV2 ^ 



l + O 



V iV2 



(2.10) 



Then, the desynchronisation bifurcation occurs for 7 = 1. We will consider values of 7 
which may be smaller than 1, but are bounded away from zero. The reason why the set 
of stationary points can be controlled in this regime is that as — > 00, the deterministic 
system x = — VV^(3;) behaves like a Ginzburg-Landau partial differential equation (PDE). 
Indeed, assume that the sites of the chain are evenly distributed on a circle of radius 
1, and that there exists a smooth function u{ip,t), G S^, interpolating the coordinates 
of x{t) in such a way that 

u(27r^,t^ = Xi{t) Vie A. (2.11) 



Then in the limit N ^ og, the discrete Laplacian in (j2.2p converges to a constant times 
the second derivative of u{-,t), and we obtain the PDE 

dtu{ip, t) = f{u{ip, t)) + ^d^^u{ip, t) . (2.12) 

Stationary solutions of (j2.12p satisfy the equation 

^u"{^) = -f{u{ip)) , (2.13) 

describing the motion of a particle of mass 7 in the inverted potential —U{ip). The 
prefactor 7 can be removed by scaling ip: Setting no(0) = u{^yJ4>), we see that uq satisfies 
the equation 

u'^ = -f{uo)=ul-uo. (2.14) 

All periodic solutions of this equation are known (cf. Section [3T3|) . and can be expressed 
in terms of Jacobi's elliptic function^ as 

^io (</>) = a(K)sn^-^==,K^ , (2.15) 

where 

• 00 is an arbitrary phase; 

• K G [0, 1) is an auxiliary parameter controlling the shape of the function: For small At, 
the function is close to a sine, while it approaches a square wave as k /" 1; 

• the amplitude a^n) is given by 

a{.f = ; (2.16) 



• the period of uo{ip) is 4vT+~k^K(k), where K denotes the complete elliptic integral 
of the first kind. 



^For the reader's convenience, we recall the definitions and main properties of Jacobi's elliptic integrals 
and functions in Appendix IXI 
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Figure 1. (a) Schematic bifurcation diagram of the hmiting PDE (|2.13p . Whenever the 
rescaled couphng intensity 7 decreases below M = 1,2, ... , a one-parameter family 

of stationary solutions bifurcates from the identically zero solution u = 0. (b) Relations 
between the shape parameter k, the amplitude a and the rescaled coupling intensity 7 for 
winding number M — 1. 



We are looking for periodic solutions of (j2.14p of period 2tt 
whenever the shape parameter k satisfies the condition 

2ttM 



4\/1 + k2K(k) 



7. Such solutions exist 



(2.17) 
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for some integer M, which plays the role of a winding number controlling the number of 
sign changes of uq. Equation (12.170 imposes a relation between shape parameter k and 
rescaled coupling intensity 7, shown in Figured] in the case M = 1. On the other hand, 
the phase 0o is completely free. 

The left-hand side of (j2.17p being bounded below by 27r, solutions of given winding 
number M exist provided 7 ^ l/M^. The smaller 7, the more different types of periodic 
solutions exist. A new one-parameter family of stationary solutions, parametrised by (po, 
bifurcates from the identically zero solution every time 7 becomes smaller than l/M^, 
M = 1,2,... (Figure IJi). 

Finally, note that for stationary points x satisfying (12. lip , with u given by ()2.15p . the 
value of the renormalised potential Vy{x)/N converges, in the limit 00, to an integral 
which can be computed explicitly (see Section 13. 3p in terms of the parameter k: 



lim 



1 



3(1 + k2 



2 + 



^ E(k) 

1 + k2 K(k) 



(2.18) 



where E denotes the complete elliptic integral of the second kind. 

If we were to add noise to the PDE (j2.12p . we would obtain a Ginzburg-Landau SPDE. 
In that case, we expect that the configurations of highest energy reached in the course of 
a typical transition from u = — lton = 1 are of the form (j2.15p . with winding number 
M = 1. As a consequence, the potential difference in ()2.18p should be governing the typical 
time of such transitions. However, proving this would involve an infinite-dimensional 
version of Wentzell-Freidlin theory, moreover in a degenerate situation, which is beyond 
the scope of the present work. We will henceforth consider situations with large, but finite 
particle number. 
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2.4 Main Results: Stationary Points for Large but Finite 

We examine now the structure of the set 5(7) for large, but finite particle number N. 
Instead of the limiting differential equation (j2.13p . stationary points satisfy the difference 
equation 

I [Xn+l - 2Xn + Xn-l] = - f {Xn) , UGA. (2.19) 

The key idea of the analysis is to interpret n as discrete time, and to consider (|2.19p as 
defining Xn+i in terms of Xn and Xn-i- Setting Vn = Xn — Xn-i allows to rewrite (I2.19P as 
the system 

Xn+l — Xn ~\~ Vn+1 5 , , 

_ (2.20) 

Vn+l = fn - 27 f{Xn) ■ 

The map (x„,fn) ^ {xn+i,Vn+i) is an area-preserving twist map ("twist" meaning that 
Xn+l is a monotonous function of Vn), for the study of which many tools are avail- 
able |Mei92j . Stationary points of the potential V-y are in one-to-one correspondence with 
periodic points of period N of this map. If we further scale u by a factor e = y^2/7, we 
obtain the equivalent map 

Xn+l =Xn + eyn+1 , 
Vn+l =yn- ej{Xn) ■ 

The regime of large particle number and finite rescaled coupling intensity 7 corresponds 
to large 7, and thus to small e. The map (|2.2ip is a discretisation of the system of 
ordinary differential equations x = y, y = —f{x), which is equivalent to the continuous 
limit equation (j2.14p . There should thus be some similarity between the orbits of the 
map (|2.2ip and of the system (j2.14p . In particular, one easily checks that the energy 

E{x,x) = ^x'^ -U{x) , (2.22) 

which is conserved in the continuous limit, changes only slightly, by an amount of order 
e^, for the map (j2.2ip (setting y = x). The map is thus close to integrable, which makes 
its analysis accessible to perturbation theory. 

Our main result, obtained by analysing the map ()2.2ip . is that the bifurcation diagram 
looks like the one shown in Figure [5J Namely, 

• For 7 > 1, and O are the only stationary points. 

• Below 7 = 1, an explicitly known number of saddles of index 1 and 2 bifurcate from 
the origin. 

• For any 2 ^ M ^ N/2, an explicitly known number of saddles of index 2M — 1 and 
2M bifurcate from the origin at 7 = 7m, where 

7M l-cos(2VjV) 1 
= — = 1 - cos(2^M/iV) = + ^ j • ^^-^^^ 

• For any fixed M, if N is sufficiently large compared to M, the above list of stationary 
points is complete for 7 > 7m- In particular, there are no secondary bifurcations of 
existing branches of stationary points, and no stationary points created by saddle-node 
bifurcation for these values of 7. 
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73 



14L 

li3(_l)i4lL30(_l)i3li4(_l)i30 
Iil(_l)i2lil(_l)illi2(_l)il 

(l^-iO(-l)^-iO)2 

I2i-l0(_l)2i-l0 
l2i(-l)2i 



— -O- 



72 



7± 



O 



/ // 



//• 

// 
// 








Figure 2. Partial bifurcation diagram for a case where = 4L is a multiple of four, 
and some associated graphs Q. Only one stationary point per orbit of the symmetry 
group G is shown. Dash-dotted curves with k dots represent fc-saddles. The symbols at 
the left indicate the zero-coupling limit of the stationary points' coordinates, for instance 
1)^-^ stands for a point whose first 2L coordinates are equal to 1, and whose last 2L 
coordinates are equal to —1. The numbers associated with the branch created at 73 are 
Li = [2L/3J , L2 = 2{L - Li), L3 = [2L/3 + 1/2J and L4, = 2(L - L3) - 1 (in case is a 
multiple of 12, there are more vanishing coordinates). 



The main difficulty is to rule out the appearance of other stationary points away from 
the origin. Indeed, for perturbed integrable maps it is easy to obtain a lower bound on the 
number of periodic points, using the Poincare-Birkhoff theorem, but it is hard to obtain 
an upper bound. One might imagine a scenario where stationary points appear far from 
the origin, which ultimately offer a more economic path for the transition from to 

We now give the precise formulation of the results. We first describe the behaviour 
between the first two bifurcation values 71 and 72. Below, gcd(a, 6) denotes the greatest 
common divisor of two integers a and 6, and Ox = {gx: g G G} denotes the group orbit 
of a point x G X under the symmetry group G. 

Theorem 2.1. There exists Ni < 00 such that when N ^ Ni and 72 < 7 < 7i = 1, the 
set S of stationary points of the potential has cardinality 

4N i3 + 2N if N is even, , 

5=3 + -, - = < 2.24 

gcdA^,2) y-i ^ 4N if N is odd , 
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There exist points A = A^j) and B = B{^) in X such that S can be decomposed a^ 

50 = Oj+ = {i+,r] , 

51 = Oa = {±A, ±RA, ±R^-^A} , 

52 = Ob = {±B, ±RB, . . . , ±R^-^B} , 

53 = Oo = {O} . (2.25) 

The potential difference between the 1-saddles and the well bottoms (which is the same for 
all 1-saddles and well bottoms) satisfies 



v,iAm-v,{i^) 1 



1 



N 4 3(1 + k2) 

where n = ^(7) is defined implicitly by 

7 



2 + ^2 ^E(k) 



4K(k)2(1 + k2) ■ 



(2.26) 



(2.27) 



The detailed proofs are given in Section [3l 

The second result, which is also proved in Section [3l concerns the behaviour for sub- 
sequent bifurcation values 7m, M ^ 2. 

Theorem 2.2. For any M ^ 2, there exists Nm < 00 such that when N ^ Nm and 
7A/+1 < 7 < 7Af , the set S of stationary points of the potential Vy has cardinality 



M 



(2.28) 



gcd(iV, 2m) ■ 

There exist points A^^^ and B^"^^ in X , m, = 1, . . . , M , such that S can be decomposed as 



So 


= Oj+ = 


'52m- 1 




S2m 




S2M+I 


= Oo = { 



m = l,...,M 
m = l,...,M 



(2.29) 



The potential difference between the saddles ^j™^ (7) and the well bottoms satisfies a similar 
relation as ()2.26p . but with k = K{m'^j). 

Remark 2.3. The proof actually yields information on the coordinates of the points 
A = ^(7) and B = B{j): 

• The coordinates of A and B satisfy the symmetries indicated in Table [H 

• If is even, the coordinates of A and B are given by 

/4K(k(7)) 



Aj{j) = 0(^(7)) sn( 
Bjil) = a(K(7))sn( 



N 



i-i),K(7)) + 



(2.30) 



where the amplitude a(fi;) is the one defined in (|2.16p . 



^If A'' is even, the orbits Oa and Ob contain A'' instead of 2A'' points, because _R^/^ 
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Figure 3. (a) Coordinates of the 1-saddles A in the case N — 32, shown for two different 
values of the coupUng 7' > 7". (b) Coordinates of the 3-saddles A^^^ in the case iV = 32, 
shown for the couphng intensities 7'/4 and 7''/4- 

• If is odd, the coordinates of A and B are given by 

^,R)=a(KR))sn(l£M))i,«R))+0(l), 

BiR) = »W7))m(i^^^^i,«(7)) + 0(1) . (2.31) 

• The components of A^'"^^ and B^^^ are given by similar expressions, with 7 replaced 
by m^7, j — ^ replaced by m{j — ^) and j replaced by mj. 

• Note that the total number of stationary points accounted for by these results is of 
the order A^^, which is much less than the 3^ points present at zero coupling. Many 
additional stationary points thus have to be created as the rescaled coupling intensity 
7 decreases sufficiently, either by pitchfork-type second-order bifurcations of already 
existing points, or by saddle-node bifurcations. However, the values 7(A^) for which 
these bifurcations occur have to satisfy limjv_^oo 7(-^) = 0. 

The existence of second-order bifurcations follows from stability arguments. For in- 
stance, for even N, the point ^(7) converges to (1, 1, ... , 1, —1, —1, . . . , —1) as 7 — > 0, 
which is a local minimum of instead of a 1-saddle. The ^-branch thus has to bifur- 
cate at least once as the coupling decreases to zero (Figured]). For odd N, by contrast, 
the point ^(7) converges to (1,1,..., 1, 0, —1, —1, . . . , —1) as 7 ^ 0, which is also a 
1-saddle. We thus expect that the point ^(7) does not undergo any bifurcations for 
^ 7 < 1 if is odd. 

2.5 Stochastic Case 

We return now to the behaviour of the system of stochastic differential equations 

d<(t) = /(xf (i)) di + I [xf+i(t) - Ix^t) + xUit)] dt + a dB,{t) , ieA. (2.32) 

Our main goal is to characterise the noise-induced transition from the configuration /~ = 
(— 1,— 1,...,— l)to the configuration = (1, 1, . . . , 1). In particular, we are interested in 
the time needed for this transition to occur, and by the shape of the critical configuration, 
i.e., the configuration of highest energy reached during the transition. 

In [BFG06ai Theorem 2.7], we obtained that in the synchronisation regime 7 > 1, 
for any initial condition xq in a ball B{I~,r) of radius r < 1/2 around any particle 
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12-L(_1)2 



12^0 



l2_Lo(-l)2-£-l ■ 



.'Aa 



yAaa 



Figure 4. Partial bifurcation diagram for a case where = 4L is a multiple of four, 
showing the expected bifurcation behaviour of the critical 1-saddlc in the zero-coupling 
limit. 



number N ^ 2 and any constant 6 > 0, the first-hitting time r+ = r'^'*(S(/+, r)) of a ball 
B{I^,r) of radius r around satisfies 

lim p-o re(iV/2-5)/.2 ^ ^{N/2+&)/a^ I ^ ^ (2.33) 

and 

lim fj^ logE^o {r+} = — . (2.34) 
(T^o 2 

This means that in the synchronisation regime, the transition between and takes a 
time of the order q^I"^" . Furthermore, for any fixed radius R G (r, 1/2), the first-hitting 
time To = t'^'*(;B(0, r)) of a ball around the origin satisfies 

lim^ P^'o {to < r+ I T+ < r_} = 1 , (2.35) 

where t_ = inf{t > t^^^^{B{I^ , R)) : xt € B{I~ , r)} is the time of first return to the small 
ball B{I~,r) after leaving the larger ball B{I~,R). This means that during a transition, 
the system is likely to pass close to the origin, i.e., the origin, being the only saddle of V-y, 
is the critical configuration of the transition. 

We can now prove a similar result in the desynchronised regime 7 < 1. 



Theorem 2.4. For 7 < 1, let 



Hi) 



y^(A(7))-n(/±) _ 1 1 



N 4 3{1 + K^) 



2 + E(k 



2- 



1 + K(/t 



+ (2-36) 



where k = ^(7) is defined implicitly by ()2.27p . Fix an initial condition xq £ B{I ,r). 
Then for any < 7 < 1, and any 6 > 0, there exists Nq{j) such that for all N > Nq{j), 



and 



lim p-o |e{2iVh(7)-5)/a2 < < ^{2Nhi^)+S)/a^ | ^ (2.37) 



lim log E^o {t+} = 2Nh{j) . (2.38) 

(7— >0 
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y-,{AW)-i/^(7±) 

N 









0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 1.3 1.4 1.5 7 



Figure 5. Value of the rescaled potential barrier height h{^) = [V-^iA) - V^[I^))/N as 
a function of the rescaled coupling intensity 7. For comparison, we also show the rescaled 
barrier height {V-y{A^'^^) — Vy{I^))/N for a stationary point of the higher winding number 
M = 2. 

Furthermore, let 

rA = r^"(U^(5Ar)) , (2.39) 

where A = A{j) satisfies ^M\f (or if N is odd). Then for any N > No{^), 

Urn ^"'{ta < t+ I r+ < t_} = 1 . (2.40) 

The relations (j2.37p and (|2.38p mean that the transition time between the synchronised 
states /~ and is of order q^^^^^)I" , while relation (12.40p implies that the set of critical 
configurations is given by the group orbit of A under G. 

The large- A'^ limit of the rescaled potential difference h{^) is shown in Figure [H The 
limiting function is increasing, with a discontinuous second-order derivative at 7 = 1. For 
small 7, h{j) grows like the square-root of 7. This is compatible with the weak-coupling 
behaviour h = (1/4 + 3/27 + 0{'y'^))/N obtained in |BFG06a] . if one takes into account 
the scaling of 7. 

The critical configuration, that is, the configuration with highest energy reached in the 
course of the transition from /~ to is any translate of the configuration ^(7) shown 
in Figure [3^. If is even, it has N/2 positive and N/2 negative coordinates, while for odd 
A^, there are (A^ — l)/2 positive, one vanishing, and (A^ — l)/2 negative coordinates. The 
sites with positive and negative coordinates are always adjacent. The potential difference 
between the 1-saddles A and the 2-saddles B is actually very small, so that transition 
paths become less localised as the particle number A^ increases, reflecting the fact that the 
system becomes translation- invariant in the large- A" limit. 
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3 Proofs 



3.1 Strategy of the Proof 

The proof of Theorems 12.11 and 12.21 is based on the fact that stationary points of the 
potential satisfy the relation 

f{Xn) + ^[xn+l-'^Xn + Xn-l] =0 , (3.1) 

where f{x) = x — x^. As mentioned in Section 12.41 this relation can be rewritten as a 
two-dimensional area-preserving twist map 

Xn+l — Xfi + Vfi+l ) 

^ ( o . ^ ) 

Vn+1 =Vn- 27~ f{Xn) ■ 

whose periodic points correspond to stationary points of the potential. In fact, we are 
going to analyse a slightly different equivalent map, which has the advantage to use the 
symmetries of the model in a more efficient way. 
The proof is organised as follows: 

• In Section [3.21 we introduce the alternative twist map, adapted to symmetries. 

• In Section [3.31 we compute the expression of the map in action-angle variables, taking 
advantage of the existence of an almost conserved quantity. 

• In Section 13.41 we compute the generating function of the map in action-angle vari- 
ables. This reduces the problem of finding periodic orbits to a variational problem 
(which is simpler than the original one). 

• The main difficulty is that the system is almost degenerate along the translation mode. 
In Section [331 we introduce Fourier variables in order to decouple the translation mode 
from the other, "oscillating" modes. 

• In Section [33l we deal with the oscillating modes, by showing with the help of Banach's 
contraction principle that for each value of the translation mode, there is exactly one 
value of the oscillating modes yielding a stationary point. 

• In Section 13.71 we deal with the translation mode, by reducing the problem to one di- 
mension, and showing that the generating function is dominated by its leading Fourier 
mode in this direction. This yields the exact number of stationary points. 

• Finally, in Section [3.81 we consider the stability of the stationary points. 



3.2 Symmetric Twist Map 



The twist map (|3.2p does not exploit the symmetries of the original system in an optimal 
way. In order to do so, it is more advantageous to introduce the variable 

Un = (3.3) 

instead of Vn- Then a short computation shows that 

Xn+l = Xn+Un-j'^fiXn) , 

(3-4) 

Un+1 = n„ - 7 [f{Xn) + f{Xn+l)\ ■ 

The map Ti : {xn,Un) i— > {xn+i,Un+i) is also an area-preserving twist map. Although 
it looks more complicated than the map ()3.2p . it has the advantage that its inverse is 
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(3.5) 



obtained by changing the sign of u, namely 

Un = Un+l + 7^"^ [f{Xn+l) + /{Xn)] ■ 

This imphes that if we introduce the involutions 

Si : {x,u) {—x,u) and 52 : {x,u) {x,—u) , (3-6) 
then the map Ti and its inverse are related by 

Ti o5i = 5i o (Ti)~^ and Ti o ^2 = ^2 o (ri)-\ (3.7) 

as a consequence of / being odd. This implies that the images of an orbit of the map 
under and ^2 are also orbits of the map. 

For large A^, it turns out to be useful to introduce the small parameter 



' ^ + 0f4,||, (3.8) 



7 V 7i7 N^fy\ \N'^ 



and the scaled variable w = u/e. This transforms the map Ti into a map T2 : {xn,Wn) 
(xn+i^Wn+i) defined by 



Wn+l =Wn- 7;^[f{Xn) + /(x„+i)] 



(3.9) 



2 

T2 is again an area-preserving twist map satisfying 

r2 o 5i = 5i o {T2)-^ and r2 o ^2 = ^2 o {T2y^ . (3.10) 

3.3 Action-Angle Variables 

For small e, we expect the orbits of this map to be close to those of the differential equation 

x = w , , ^ 

w = -f{x) , 

which is equivalent to the second-order equation x = —f{x) describing the motion of a 
particle in the inverted double- well potential —U{x), compare (|2.14p . Solutions of (13.1ip 
can be expressed in terms of Jacobi elliptic functions. Indeed, the function 

C{x,w) = ^{x^ + w'^)-^x'^ (3.12) 
being a constant of motion, one sees that w satisfies 



w = ±VHC? - x2)(6(C7)2 - x2)/2 , (3.13) 

where 

a(C)2 = 1 - Vl - 4C , 

b{Cf = 1 + VI -4C . (3.14) 



14 



This can be used to integrate the equation x = w, yielding 



bjC) 
V2 



t = F I Arcsin 



V«(c) 



,k{C)] , 



(3.15) 



where k(C) = a{C)/b{C), and F(0, k) denotes the incomplete elliptic integral of the first 
kind. The solution of the ODE can be written in terms of standard elliptic functions as 



x{t) = a{C)sn(^^^t,K{C) ] , 



w{t) 



/2Ccn 



(^t,.(C))dn 



(3.16) 



We return now to the map T2 defined in (j3.9p . The explicit solution of the continuous- 
time equation motivates the change of variables : {x, w) 1— > {(p, C) given by 



F (^Arcsin ( 



h{c) V 



\a{C) 



,k{C)] , 



C = \{x' + w')-\x\ 



(3.17) 



One checks that <I>i is again area-preserving. The inverse is given by 
x = a(C)sn(^V9,K(C7)) , 



w 



'2Ccn 



\ V2 



(3.18) 



The elliptic functions sn, cn and dn being periodic in their first argument, with period 
4K(k), it is convenient to carry out a further area-preserving change of variables $2 : 

((/?,C) ^ (ip,!), defined by 



where 



nic) 



V' = n{c)ip 

b{C) 



I = h{C) , 



TT 



h{C) 



c 



^) 



(3.19) 



(3.20) 



V2 2K(k(C)) ' 

Using the facts that C and h = h(C) can be expressed as functions of k = n{C) by 
C = K^jiX + ^2)2 and 6^ = 2/(1 + k^), one can check that 



h{C) 



4 (1 + k2)e(k)-(1-k2)k(k) 



37r 



(1 + ^2)3/2 



k=k(C) 



0, 



2\/2' 
"3^ 



(3.21) 



We denote hy ^ = ^2° ^\ the transformation (x, vS) 1— > (V', and hy T = ^ 0T20 ^ ^ 
the resulting map. 



Proposition 3.1. The map T = T[e) has the form 

V-n+l = V'n + en{In) + £^f{lpn, In, e) 
In+1 = In + £^g{tpn, In,£) , 



(mod 2tt) , 



(3.22) 
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where = Q{h ^{I))- The functions f and g are Tr-periodic in their first argument, 

and are real-analytic for ^ / ^ /i(l/4) — 0{£^). Furthermore, T satisfies the symmetries 



T o Si = Si o T ^ and T o S2 = S2 o T ^ , 
where Si(^/;, /) = {—tp, I) and T,2{ip, /) = (vr — ip, I). 



(3.23) 



Proof: First observe that $1 and ^ are analytic whenever (x,w) is such that C < 1/4. 
The map T will thus be analytic whenever {ipn,In) is such that C{xn,Wn) < 1/4 and 
C{xn+i,Wn+i) < 1/4. A direct computation shows that 



C{Xn+l,Wn+l) — C{x 



XnWn + 2Xnw1^ - 'ix^Wn + ^X^Wn 



(3.24) 



This implies that I^+i = In + 0(£^), and allows to determine g{ip, /, 0). It also shows that 
T is analytic for /„ < /i(l/4) — 0{£^). Furthermore, writing a„ = a{C{xn,Wn)), we see 
that (|3.24|) implies a^+i — an = 0{£^) and similarly for 6„, k„. This yields 



ip{Xn+l,Wn+l) - ip{x 



du 



bn Jx„/a„ y/(1 - K2n2)(l - U^) 



+ 0(6=^) 



(3.25) 



which implies the expression for ipn+i- We remark that the fact that T is area-preserving 
implies the relation 



1 



d{lpn+lJn+l) 
d{tpn,In) 



1 + 



d^f{^J,Q)+dig{iP,I,0) 



+ 0(e^) 



(3.26) 



which allows to determine f{ip,I,Q). The fact that / and g are vr-periodic in their first 
argument is a consequence of the fact that T2{—x,—w) = —T2{x,w). Finally the rela- 
tions (j3.23p follow from the symmetries (j3.10p . with Sj = <I> o 5j o $^1. □ 



A perturbation expansion at / = shows in particular that 

m = i-7i+o{i'). 



(3.27) 



An important observation is that i7(/) is a monotonously decreasing function, taking 
values in [0, 1]. The monotonicity of Q makes T a twist map for sufficiently small e, which 
has several important consequences on existence of periodic orbits. 

We call rotation number of a periodic orbit of period N the quantity 



1 



27riV 



N 



^(V'n+1 - tpn) (mod 2tt) 



■n=l 



(3.28) 



Note that because of periodicity, is necessarily a rational number of the form u = M/N, 
for some positive integer M. We denote by the set of points ip in the torus 
satisfying (j3.28p . It is sometimes more convenient to visualise as the set of real 
A''-tuples {ipi, . . . , iPn) such that 



ipi < ip2 < • • • < i'N < "01 + 27rA^i/ . 



(3.29) 
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In the sequel, we shall use the shorthand stationary point with rotation number u instead 
of stationary point corresponding to a periodic orbit of rotation number v. 
The expression p.22p for T implies that 

0(/o) 



27r 



(3.30) 



The following properties follow from the Poincare-Birkhoff theorem, whenever e > is 
sufficiently small: 

• For each positive integer M satisfying 



(3.31) 



the twist map T admits at least two periodic orbits of period N and rotation number 
V = M/N . Note that Condition (|3.31|) is compatible with the fact that O bifurcates 
for7 = 7M, M = l,2,...,LiV/2j. 

Any periodic orbit of period N of the map T is of the form 

V'o + 27ri/n + 0(e2) , 



(3.32) 



for some ipQ and some v = M/N, where M is a positive integer satisfying ()3.3ip . 
Going back to original variables, we see that these periodic orbits are of the form 



On sn 



y2Cn 



cn 



2K(k„) 



vr 



V'n, K-n dn 



(3.33) 



where a„ = a{Cn), k„ = n{Cn) and 



Cn = +0{e)= ( M Vt ) + 0{e) 



(3.34) 



This allows in particular to compute the value of the potential at the corresponding sta- 
tionary point. 

Proposition 3.2. Let e > be sufficiently small, and let x* be a stationary point of the 
potential V-y, corresponding to an orbit with rotation number u = M/N . Then 

V^ix") _ 1 



2 + _ ^^{k) 

1 + ^2 



K(k) 



N 3(1 + k2) 

where k = k{C), and C satisfies Q.{C)'^ = M^^. 
Proof: The expression (j2.4p for the potential implies that 

N N 

- ^([/(x„) + -wl + 0{e^)) = - Y,{< -Cn + 0{e')) 



(3.35) 



N 



n=l 

N 



-y 

N ^ 



n=l 



c 



2cn^ 
cn^ 







2 

2K{k) 
vr 

2K{k) 
vr 



V'n, K-j dn 
K I dn' 



N 

n=l 

2K{k) 

IT 

2 1 2K(k) 



l + 0(e] 



vr 



(3.36) 
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where C = ^(My^) and k = k(C). The integral can then be computed using the 
change of variables 2K(K)'(/;/7r = F{(f), k). Finally, recall that C = k? /{I + k?)"^ . □ 



One can check that V^{x*)/N is a decreasing function of k, which is itself a decreasing 
function of M"^^. As a consequence, V-f{x*)/N is increasing in M'^j. This implies in 
particular that the potential is larger for larger winding numbers M. 

Remark 3.3. The leading term in the expression (13.350 for the value of the potential 
is the same for all orbits of a given rotation number v. Since stationary points of the 
potential of different index cannot be at exactly the same height, the difference has to 
be hidden in the error terms. In |BFG06a] . we showed that near the desynchronisation 
bifurcation, the potential difference between 1-saddles and 2-saddles is of order (1 — 7)^/^. 
For large A^, we expect this difference to be exponentially small in 1/N, owing to the 
fact that near-integrable maps of a form similar to (j3.22p are known to admit adiabatic 
invariants to that order (cf. [BK961 Theorem 2]). 



3.4 Generating Function 

In this section, we transform the problem of finding periodic orbits of the near-integrable 
map T into a variational problem. The fact that T is a twist map allows us to express /„ 
(and thus In+i) as a function of ipn and ipn+i- A generating function of T is a function 
G{ipn,''Pn+i) such that 

diG{lpn,tpn+l) = -In , d2G{lp„, 1pn+l) = In+1 ■ (3.37) 

It is known that any area-preserving twist map admits a generating function, unique up 
to an additive constant. Since T depends on the parameter e, the generating function 
G naturally also depends on e. However, we will indicate this dependence only when we 
want to emphasise it. 



Proposition 3.4. The map T admits a generating function of the form 



G(V'1,V2) =£^0 



-'^) +2e=^f;Gp(^^^^^^,e)cos(p(V'i + ^2)) , (3.38) 



where the functions GQ{u,e) and Gp{u,e) are real- analytic for u > ©(l/jlog e|), and satisfy 

G',iu,o) = n-\u) , 

Gp{u,0) = — J g{'il^,n (n),0) sin(-2pV) dV' • (3.39) 

Proof: Fix (^^2, ^2) = T{^pi,Ii). The fact that T{ipi + vr, h) = (^^2 + vr, h) implies 

G(V'i+vr,V'2 + vr) = G(V'i,V'2) + c (3.40) 
for some constant c. If we set G(^i, -02) = G{ip2 — ipij'ipi + tp2), we thus have 

G{u,v + 2tt) =Giu,v)+c. (3.41) 

This allows us to expand G as a Fourier series 

00 

G(Vi,V2)= V Gp(^2-V'i,e)e'^(^^+^^)+;f (Vi+^2) . (3.42) 

p= — 00 
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Next we note that the symmetry (|3.'23p implies T{—ip2,l2) = i—tpi,h), and thus 

diG{i>u^2) = -d2G{-^P2,-^Pi) . (3.43) 

Plugging (j3.42p into this relation yields 

c = and G-p{u,e) = Gp{u,e) , (3.44) 

which allows to represent G as a real Fourier series as well. Computing the derivatives 
h = -9iG('0i,'02) and I2 = 52G(V'i,V'2) yields 

00 

I2-Ii=2 ipGp(V2-^i,e)e^^('^^+'^^) , (3.45) 

p=— 00 

which shows in particular that Gp{u,e) = 0{e^) for p ^ 0, as a consequence of (j3.22p . 
This implies h = Go(V'2 — '01, e) + C>{£^), and thus u = GQ{e^}{u) + C'(e^),e). Renaming 
Go{u,s) = sGo{u/£,£) and Gp{u,e) = £^Gp{u/e,e) yields (|3.38p . Evaluating (|3.45p for 
e = and taking the Fourier transform yields the expression (j3.39p for Gp{u,0). □ 

The relations (j3.39p allow to determine the expression for the generating function of 
the map T, given by (I3.22p . In particular, one finds 

Go(u,0) = uTi~\u) - n~^{u) , (3.46) 

so that 

Go{n{C),0) = h{G)n{G) - C 



1 



3(1 + ^2 



2 + ^E(k) 

1 + ^2 K{k) 



(3.47) 



with K = k(G). Note that this quantity is identical with the leading term in the ex- 
pression (j3.35p for the average potential per site. This indicates that we have chosen the 
integration constant in the generating function in such a way that V-y and Gn take the 
same value on corresponding stationary points. 

The main use of the generating function lies in the following fact. Consider the A^-point 
function 

GNiiJi, ...,M = G(Vi, V2) + ^(^^2, V's) + • • • + G{^N,i^i + 2t:Nv) , (3.48) 

defined on (a subset of) the set T^. The defining property ()3.37p of the generating function 
implies that for any periodic orbit of period N of the map T, one has 

d 



GN{i^i,...,i^N) = -In + In = 0. fov n = I, . . . , N . (3.49) 



In other words, A^-periodic orbits of T with rotation number 1/ are in one-to-one corre- 
spondence with stationary points of the A^-point function Gn on T^. 

The symmetries of the original potential imply that the A^-point generating function 
satisfies the following relations on : 

Gn{iPi, ■ ■ ■,'4'n) = Gn{iP2, ■ ■ .,ipN,'4'i + 2ttNi^) , 
GNiipi, ■ ■ ■ An) = GNi-ipN, . . . , -ipi) , 

Gn{iPi, ■ ■ ■,i'N) = Gn{iPi + vr, . . . ,'0Af + tt) . (3.50) 
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At this point, we are in the following situation. We have first transformed the initial 
problem of finding the stationary points of the potential into the problem of finding 
periodic orbits of the map Ti , or, equivalently, of the map T. This problem in turn has 
been transformed into the problem of finding the stationary points of Gat. Obviously, 
the whole procedure is of interest only if the stationary points of Gjy are easier to find 
and analyse than those of V^. This, however, is the case here since the A^-point function 
is a small perturbation of a function depending only on the differences V'n+i ~ V'n- In 
other words, Gn can be interpreted as the energy of a chain of particles with a uniform 
nearest-neighbour interaction, put in a weak external periodic potential. 



3.5 Fourier Representation of the Generating Function 

The main difficulty in analysing the stationary points of the A^-point generating function 
G AT comes from the fact that it is almost degenerate under translations of the form ^„ i— > 
V'n + c Vn. The purpose of this section is to decouple the translation mode from the other 
variables, by introducing Fourier variables. 

We fix = M/N . Any stationary point of Gn on admits a Fourier expansion of 
the form 

N-l 

ll^n = 27TUn + V'g^^^" , (3.51) 

q=0 



where cu = e'^'^ ' , and the Fourier coefficients are uniquely determined by 

N 

N 



I ^ 

^9 = ^1^ ^"''"(V'n - 2TTun) = 7P_g . (3.52) 

n=l 

Note that ipg = V'g+TV for all q. Stationary points of Gj\f correspond to stationary points of 
the function Gn, obtained by expressing Gn in terms of Fourier variables (ipo, • • • , V'A^-i)- 
In order to do so, it is convenient to write 

= A + £ an{^pl, . . .,Vn-i) , 

i^n + V'n+i = 24'o + 2TTiy{2n + 1) + e^Pnii'i, ■ ■ ■,i^N-i) , (3.53) 
where A = 2'kv je and 

a. 



^ 9=1 

N-l 



9=1 

Note that On is of order 1 in e for any stationary point because of the expression (|3.22|) 
of the twist map. Taking the inverse Fourier transform shows that \tpq{u!'^ — 1)| = 0{e'^) 
and iV'gl = O(e^) for g 7^ 0, and thus /3„ is also of order 1. 
Expressing Gat in Fourier variables yields the function 

oo 

GN{4^o,---,i'N-i)= ""^'"^^ 9p{k^---AN-i) , (3.55) 

p=— oo 
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where (we drop the e-dependence of Gq and Gp) 

N 

goiipi, . . .,ipN-i) = Go(A + e^a„) , 

n=l 
N 

gp{4>i, i^N-i) =e'Yl + e'«n)^^*'('"+') e'^'P"- for p / . (3.56) 

n=l 

We now examine the symmetry properties of the Fourier coefficients Qp. Table [2] shows 
how the Fourier variables transform under some symmetry transformations, where we only 
consider transformations leaving invariant. As a consequence, the first two symmetries 
in (j3.50p translate into 

Qpi^l^i, . . . ,i)N-i) = u?^^'^ gp{uj'^i, . . . ,uj'^~'^ipN-i) , 

gp{ipi, . . . ,ipN-i) = u]''^'^'^'^ g-p{-uj^~^^N-i, • • • , -t^V'i) ■ (3.57) 

We now introduce new variables Xg, g 7^ 0, defined by 

X, = -iu;-'''^°/2">, = -X^. (3.58) 

The Xq S'i'G defined in such a way that they are real for stationary points satisfying, in 
original variables, the symmetry Xj = —Xn^-j for some uq. For later convenience, we 
prefer to consider q as belonging to 

a-{- 

rather than {1, . . . , — 1}. We set x = {Xq}qeQ a-iid 

00 

= E e''^'^»5p(^o,x), (3.60) 

p=—oo 

where 

5p(^o,x) = gpiii'q = i^^''^''/'""xages) • (3.61) 
Lemma 3.5. The function G]\f{iJjQ,x) 2tiv -periodic in its first argument. 
Proof: By (|3.57p . we have 

~gp{^o + 2vrz., x) = uj'^^^' U^o, x) (3.62) 
Since q^^P''^'^" = cj^^^, replacing TpQ by ■(/'o + 27rz^ in (j3.60p leaves Gat invariant. □ 



R 


Xj 1 ^ Xj-\-\ 




V'n+l 




-> Uj'^i)q + 2TTu6gO 


CS 


Xj ^ -XN+l-j 




—1pN+l-n 


i>q ^ 


-UJ-'i^N-q - 2'KV{N + l)(5gO 


c 


Xj 1 ^ Xj 






i'q ^ 





Table 2. Effect of some symmetries on original variables, angle variables, and Fourier variables. 



A^- 1 




A 


2 


, . . . , 





\{o} 



(3.59) 
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Since also has period vr, it has in fact period 

, K = gcd{N,2M) . (3.63) 

Our strategy now proceeds as fohows: 

1. Show that for each ipQ, and sufficiently small e, the equations dGj^/dxq = 0, g G Q, 
admit exactly one solution x = X*(V'o)- This is done in Section [3.61 with the help of 
Banach's fixed-point theorem. 

2. Show that for x = X*{'^o)y the equation dG]\f/dipo = is satisfied by exactly ^N/K 
values of ^o- This is done in Section 13.71 by estimating the Fourier coefficients of 
BGn /di^o with the help of complex analysis. 



3.6 Uniqueness of x 

In this section, we show that the equations 

dGN 



dXq 



, q€Q, (3.64) 



admit exactly one solution x = X*(V'o) for each value of ipQ. The proof is based on a 
standard fixed-point argument: First we show in Lemma 13.61 that (j3.64p is equivalent 
to the fixed-point equation p = Tp for a quantity p related to x- Then we show in 
Proposition 13.81 that T is contracting in an appropriate norm, provided £ is sufficiently 
small. 

It is useful to introduce the scaled variables 

2 

Pq = Pgix) = -^Xq sm{TTq/N) (3.65) 
and the function T^^''^\p), p = {pq}qeQ, defined for £ G Z and a, 6 ^ by 

a b 

r'f''\p)= E 1{E, q',=i} n n tan(vrg7iV) ^ " ^^'^^^ 

q'l,-,q',eQ 

By convention, any term in the sum for which q'j = N/2 for some j is zero, that is, we set 
1/ tan(7r/2) = 0. A few elementary properties following immediately from this definition 
are: 

• rf '°V) = Sio; 

• rf'''\p) = for \i\ > (a + b)N/2; 

• Upq = for q^KZ, then v'f'''^ (p) =Oioi £^ KZ; 

. If = p.q for all q, then Tf''\p') = (-l)TLt^(p). 

The following result states that the conditions p.64p are equivalent to a fixed-point 
equation. 



Lemma 3.6. Let 
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with the convention that -f/pjv/2(^) = ^'^(A). Then the stationarity conditions p.64p are 
fulfilled if and only if p = p{x) satisfies the fixed-point equation 

p = Tp = p^''^ + ^p,e) , (3.68) 

where the leading term is given by 
1 



4°^ 



ifq^KZ, 



and the remainder is given by ^q{p,£) = ^'^g\p,e) + ^^^\p,e), with 

T . ^2(a+fe) , , 

(3.70) 

The proof is a straightforward but lengthy computation, which we postpone to Ap- 
pendix [B1 

Note the following symmetries, which follow directly from the definition of p^'^^ and 
the properties of F^"' : 

• For all q £ Q, = Pq^\ because -ff_p^_g(A) = Hp^g{A), and thus pq^^ E M; 

• If pq = OfoT q^ KZ, then ^q{p, e) = for g JTz'; 

• If = p-q for all q, then $g(p',e) = $_g(p,e); 

Remark 3.7. The condition kN + g G 2MZ , appearing in the definition of p^^\ can only 
be fulfilled if g G iVZ + 2MZ = KZ . If this is the case, set = nK, 2M = mK, q = iK, 
with n and m coprime. Then the condition becomes mp — kn = i. By Bezout's theorem, 
the general solution is given in terms of any particular solution (pQ, ko) by 

p = Pq + nt , k = ko + mt , t £ Z . (3.71) 

Thus there will be exactly one p with 2\p\ < n. If A is very large, and M is fixed, then 
n = N/K is also very large. Since the Gp(A), being Fourier coefficients of an analytic 
function, decrease exponentially fast in the sum in (j3.69p will be dominated by the 
term with the lowest possible \p\. 

We now introduce the following weighted norm on C ^: 

||p||A = supe^l«l/2A-^|p,| , (3.72) 

where A > is a free parameter. One checks that the functions Go(^) and Gp(A) are 
analytic for ReA > 0(1/ log|e|). Thus it follows from Cauchy's theorem that there exist 
positive constants Lq, r < A — ©(l/loglej) and Aq such that 



\G^^\^)\ ^ Lo^ and |G('^)(A)| ^ Lo^e-^°H (3.73) 
for all a ^ and p G Z . For sufficiently small e, it is possible to choose r = A/2. 
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Proposition 3.8. There exists a numerical constant ci > such that for any A > 0, any 
A < Ao and any Rq > ciLo[A|Gg(A)|]~-'^, there is a strictly positive Eq = eo(A, A, Aq, -Ro) 
such that for all e < Eq, the map T admits a unique fixed point p* in the ball 13x{0, Rq) = 
{p € C^: IIpIIa < Rq}- Furthermore, the fixed point satisfies 

• p* = whenever q KT, ; 

• p*_q = p*, and thus p* gM. for all q. 

The proof is again a straightforward but lengthy computation, so we postpone it to 
Appendix [BJ 

A direct consequence of this result is that for any ipQ, and sufficiently small e, there is 
a unique p* = /5*(^o) (and thus a unique x*{'4'o)) satisfying the equations dG]\f/dxq = 
for all q £ Q. Indeed, we take i?o sufficiently large that our a priori estimates on the Xq 
imply that p £ Bo{0,Ro). Then it follows that p is unique. Furthermore, for any A < Aq, 
making e sufficiently small we obtain an estimate on ||p*||a- 

3.7 Stationary Values of ^/^o 

We now consider the condition dGN/dtpo = 0. As pointed out at the end of Section [33} 
Giv('0Oi x) is a vrK/A^-periodic function of ipQ. For the same reasons, x*(V'o) is also ttK/N- 
periodic. Hence it follows that the function tpQ ^ Gn{'^q-,X*{'4'o)) has the same period as 
well, and thus admits a Fourier series of the form 

oo 

GN{^o,x''m)= E fffce^i^'^oA^/^, (3.74) 

k=—oo 

with Fourier coefficients 

E ^''^"-'''"'^^'U^o^xim'^i^^ (3.75) 



9k 



p=—oo 



(we have chosen [0, 27r] as interval of integration for later convenience). Using the change 
of variables ipQ —ipo in the integral, and the various symmetries of the coefficients (in 
particular (j3.57p ). one checks that g^k = g^- Therefore (j3.74p can be rewritten in real 
form as 

oo 

Giv(V'o,X*(^o))=5o + 2E 

gkCos(2kTpQN/K) . (3.76) 

k=l 

Now BGn /diljQ vanishes if and only if the total derivative of Gn{iPo, X*(V^o)) with respect to 
■00 is equal to zero. This function obviously vanishes for ipo = lirK/lN, i = 1, . . . , AN/K, 
and we have to show that these are the only roots. 

We first observe that the Fourier coefficients g^ can be expressed directly in terms of 
the generating function (13.380 . written in the form 

G{u,v,e) = eGo{u,e) + E^Y.^pi'^^^)'^'^'^'' ■ (3-77) 
In the sequel, a^{tpo) and /3^(V'o) denote the quantities introduced in (I3.54p . evaluated at 
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Lemma 3.9. The Fourier coefficients are given in terms of the generating function by 



N 

9k = 

where 



9k = i^ r e-^i'^^o^/^ Ao(^o) dV'o , (3.78) 



Ao(V^o) = G(A + e2a*(^o),V;o + 7ri/ + ie2^5(V'o),e) • (3.79) 
Proof: The coefficient g^ can be rewritten as 

1 ^ r2TT 

where 

+ e3^pe2iP^«Gp(A + e2a^(^o))^pA^{2n+i)ei^'p/3n(^o) . (3,31) 

Using the periodicity of x*, one finds that a* (V'o + ^iriy) = a*^;^('(/'o) and similarly for /?*, 
which implies An{ipo) = Ao(V'o + 27rz^n). Inserting this into p.80p and using the change 
of variables ipo ipo — 2'ku in the nth summand allows to express g^ as the {2kN / K)t\]. 
Fourier coefficient of Aq. Finally, Ao(V'o) can also be written in the form (I3.79p . □ 

Relation (j3.78p implies that the cjk decrease exponentially fast with A:, like Q-'^^okN/K ^ 
Hence the Fourier series (I3.76P is dominated by the first two terms, provided N is large 
enough. In order to obtain the existence of exactly AN/K stationary points, it is thus 
sufficient to prove that ^1 is also bounded below by a quantity of order q~'^'^oN/k ^ 

Proposition 3.10. For any A > 0, there exists £i(A) > such that whenever e < ei(A), 

sign(5i) = (-l)i+2M/i. _ (3,82) 

Furthermore, 

^^exp|-^^Ao(A)^l yk^2, (3.83) 



\gi\ ^ 4 'K 

where Ao(A) is a monotonously increasing function of A, satisfying Ao(A) = ^/2tt/S. + 
O(A^) as A \ 0, and diverging logarithmically as A 1. 

Proof: First recall that A = 2-kv je = 2'KM/Ne, where M is fixed. Thus taking e smah 
for given A automatically yields a large N. Combining the expression (I3.22P for the 
twist map and the defining property (j3.37p of the generating function with the relations 
u = (tpn+i — V'n)/e and v = i/jn + V'n+1; one obtains the relation 

3 

dyG{u,v,e) = ^[g(^^{v - eu),n~\u),e) + 0{e^)\ . (3.84) 
It follows from ([X^ and the definition ([X^ of h{C) that 

5(V', /, 0) = ^ [1 + 2«;2 _ + 3x'] 

^ '^"'■[l + 4C-6x2 + 4x^] , (3.85) 



n{I) 4 
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Figure 6. The integration contour F used in the integral p.89p . 



where x and w have to be expressed as functions of -0 and / via ()3.17p and ()3.19p . In 
particular, we note that 



/2C n 



^ _ of-n^ 

a 2K(K)dV ^ Ui; 



where we used (|3.20p again. This allows us to write 

1 d 

8d0 



3x^ + 



(3.86) 



(3.87) 



A similar argument would also allow to express the first-order term in e of g{il^, I, e) as a 
function of x = x{'4j,I). Also note the equality 



A = n{I) + 0{e) 



TrbjC) 
2V2K(k] 



+ Oie) 



1 



+ Oie) , 



(3.88) 



which follows from the relation (j3.30p between and 0(7). 

The properties of elliptic functions imply that for fixed I, ip ^ x^ip, I) is periodic in 
the imaginary direction, with period 2Ao = 7rK(\/l — k^)/K(k), and has poles located 
in -0 = nvr + (2m + l)iAo, n,m G Z. As a consequence, the definition of the map 
T = ^> o T2 o <I>~^ implies in particular that g{ip, I, e) is a meromorphic function of ip, 
with poles at the same location, and satisfying g{Tp + 2iAo,/, e) = g{tp,I,e). These 
properties yield informations on periodicity and location of poles for Ao(V'o)) in particular 

Ao(Vio + 2iAo) =Ao(V^o) + 0(e')- 

Let r be a rectangular contour with vertices in — vr/2, —tt/2 — 2i Aq, 37r/2 — 2 i Aq, and 
37r/2, followed in the anticlockwise direction (Figured]), and consider the contour integral 



J 



1 

2^ 



-2ikzN/K ^^^^^ dz. 



(3.89) 



The contributions of the integrals along the vertical sides of the rectangle cancel by peri- 
odicity. Therefore, by Lemma [3.91 and the approximate periodicity of Aq in the imaginary 
direction. 



J 



1 

'n 



gk 



-2kXoN/K 



gk + OiNe' 



On the other hand, the residue theorem yields 



J = 27r i ^ g~2i/c.,iv/i^ Res(Ao(z), zj) , 



(3.90) 



(3.91) 
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where the Zj denote the poles of the function Ao(2;), lying inside T. There are two such 
poles, located in zi = — i Aq— vri/+eA+C'(e^), and Z2 = zi+tt, and they both yield the same 
contribution, of order Q-^okN(i+o{e ))/K ^ ^j^g sum. Comparing (j3.90p and (j3.9ip shows 
that Qk/N is of the same order. Finally, the leading term of can be determined explicitly 
using ()3.87p and Jacobi's expression (|A.lip for the Fourier coefficients of powers of elliptic 
functions, and is found to have sign (—1)1+2*^^ for sufficiently large N . Choosing e small 
enough (for fixed eN) guarantees that gi dominates all for k ^ 2. □ 

Corollary 3.11. For e < ei{A), the N -point generating function Gn admits exactly 
AN/K stationary points, given by -00 = ^ttK/2N , i = 1, . . . , AN/K , and x = X^H^o)- 

Proof: In the points ipo = £ttK/2N ^ the derivative of the function i/'o ^ Giyiipo, X*(V'o)) 
vanishes, while its second derivative is bounded away from zero, as a consequence of 
Estimate ()3.83p . Thus these points are simple roots of the first derivative, which is bounded 
away from zero everywhere else. □ 



3.8 Index of the Stationary Points 

We finally examine the stability type of the various stationary points, by first determin- 
ing their index as stationary points of the A^-point generating function G^, and then 
examining how this translates into their index as stationary points of the potential V^. 

Proposition 3.12. Let (V'O; X*(V'o)) be a stationary point of Gn with rotation number 
V = N/M . Let X* = x*{'iPq) be the corresponding stationary point of the potential V^, and 
let K = gcd(iV, 2M). 

• If2M/K is odd, then the points x*(0), x*{Ktt/N), . . . are saddles of even index ofV^, 
while the points x'^{Ktt/2N), x*{3Ktt/2N), . . . are saddles of odd index ofV^. 

• If2M/K is even, then the points 2;*(0), x*{Ktt/N), . . . are saddles of odd index ofV^, 
while the points x*(K'it/2N), x*{3Kit/2N), . . . are saddles of even index ofV-y. 

Proof: We first determine the index of {tpo, x*(V'o)) as stationary point of Gat. Using the 
fact that Gq(A) is neg ative {n (A) being decreasing), one sees that the Hessian matrix 
oi Gn is a small perturbation of a diagonal matrix with — 1 negative eigenvalues. The 
A^th eigenvalue, which corresponds to translations of tpo, has the same sign as the second 
derivative of ipo GAr(V'o, x*(^o)), which is equal to (— 1)^*^/^ signcos(2';/;oA^/-fC). Thus 
{ipo, x*(V'o)) is an A^-saddle of Gn if this sign is negative, and an (A^ — l)-saddle otherwise. 
The same is true for the index oi = {ipi, . . . , ■i/'Af) as a stationary point of Gat. 

Let R be the so-called residue of the periodic orbit of T associated with the stationary 
point. This residue is equal to (2 — Tr(L'T^))/4, where DT^ is the Jacobian of at 
the orbit, and indicates the stability type of the periodic orbit: The orbit is hyperbolic if 
-R < 0, elliptic if < i? < 1, and inverse hyperbohc if i? > 1. It is known |MM83| that the 
residue R is related to the index of ijj be the identity 

ii = -l^i!5£=^^ML. (3.92) 

4n7.i(-9i2G(*„vi+i)) 

In our case, —di2G{tpj,tpj^i) is always negative, so that R is positive if ip is an (A^ — 1)- 
saddle, and negative if ip is an A^-saddle. 

Now x*{'iPq) also corresponds to a periodic orbit of the map (j3.2p . whose generating 
function is Xn+i) = — x„+i)^ + ^f/(x„). The corresponding A^-point generating 



27 



function is precisely (2/7) V^. Since the residue is invariant under area-preserving changes 
of variables, we also have 



1 det(Hessy^(x*)) 
27nf=i(-5i2^(x*,x^Vi)) 



^ = -t^ ^/7T ::T . ■ (3.93) 



In this case, the denominator is positive. Therefore, HessV^(x*) has an even number of 
positive eigenvalues if ip is an A^-saddle, and an odd number of positive eigenvalues if ■0 is 
an {N - l)-saddle. □ 

We can now complete the proofs of Theorem 12.11 and Theorem 12.21 

Proof of Theorem 12. 1[ We first recah the following facts, established in |BFG06aj . 
Whenever 7 crosses a bifurcation value ■jm, say from larger to smaller values, the index 
of the origin changes from 2M — 1 to 2M + 1. Thus the bifurcation involves a centre 
manifold of dimension 2, with 2M— 1 unstable and N —2M —1 stable directions transversal 
to the manifold. Within the centre manifold, the origin repels nearby trajectories, and 
attracts trajectories starting sufficiently far away. Therefore, all stationary points lying 
in the centre manifold, except the origin, are either sinks or saddles for the reduced two- 
dimensional dynamics. For the full dynamics, they are thus saddles of index (2M — 1) or 
2M (c.f. |BF(Tn6a[ Section 4.3]), at least for 7 close to jm- 

We now return to the twist map in action-angle variables (j3.22p . The frequency 0,(1) 
being maximal for / = 0, as e increases, new orbits appear on the line 1 = 0, which 
corresponds to the origin in x-coordinates. Orbits of rotation number u = M/N can only 
exist if £0(0) = e ^ 2ttv + O(e^), which is compatible with the condition 7 < 7^/ . 

Consider the case of a winding number M = 1, that is, of orbits with rotation number 
u = 1/N, which are the only orbits existing for 72 < 7 < 71- We note that 7 > 72 
implies A = 271 /Ne > 1/2 — 0{e), and thus there exists A^i < 00 such that the condition 
N ^ Ni automatically implies that e is small enough for Corollary 13.111 to hold. Now, 
Proposition 13.121 yields: 

• If is even, then K = gcd(A^, 2) = 2, and there are 2N stationary points. The points 
x*(0),a;'^(27r/A^), . . . must be 2-saddles, while the points x*{it/N),x*{3tt/N), . . . are 

1- saddles; 

• If is odd, then K = gcd{N, 2) = 1, and there are AN stationary points. The points 
x*{0),x*{tt/N), . . . must be 1-saddles, while the points x*{7r/2N), x* {3tt/2N), ... are 

2- saddles. 

Going back to original variables, we obtain the expressions (I2.30p and (|2.3ip for the co- 
ordinates of these stationary points. The fact that they keep the same index as 7 moves 



away from 74/ is a consequence of Relation (I3.93P and the fact that the corresponding sta- 
tionary points of Gat also keep the same index. Finally, Relation (j2.26p on the potential 
difference is a consequence of Proposition 13. 2[ This proves Theorem 12.11 □ 

Proof of Theorem 12. 2[ For larger winding number M, one can proceed in an analo- 
gous way, provided is sufficiently large, as a function of M, for the conditions on e to 
hold. This proves Theorem 12.21 □ 

Finally, Theorem 12.41 is proved in an analogous way as Theorems 2.7 and 2.8 in 
[BFGn6a| . using results from |FW98| (see also [KifSH [Sug96l ) . 
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A Jacobi's Elliptic Integrals and Functions 

Fix some k G [0,1]. The incomplete elliptic integrals of the first and second kind are 
defined, respectively, b}jl 

f'^ dt I 

F((/),k)= / —^^^^=, E((/),k)= / V 1 - k2 sin^ t dt . (A.l) 

h v 1 - k2 sin^ t Jo 

The complete elliptic integrals of the first and second kind are given by 

K(k) = F(^/2,k) , E(k) = E(7r/2,K) . (A.2) 

Special values include K(0) = E(0) = tt/2 and E(l) = 1. The integral of the first kind 
K(k) diverges logarithmically as k 1- 

The Jacobi amplitude am(u, k) is the inverse function of F(-, k), i.e., 

(/) = am(n, k) 44> u = F((j),K). (A. 3) 

The three standard Jacobi elliptic functions are then defined as 

sn(n, k) = sin(am(u,K)) , 

cn(n, k) = cos(am(ti, k)) , (A. 4) 



dn(n, k) = \J 1 — K? sin (am (u, k)) . 
Their derivatives are given by 

sn'(M, k) = cn(ti, k) dn(ti, k) , 

cn'(n, k) = — sn(ii, k) dn(n, k) , (A.5) 
dn'(u, k) = —h? sn(n, k) cn(n, k) . 

The function sn satisfies the periodicity relations 

sn(ti + 4K(k), k) = sn(n, k) , 
sn(n + 2 i K(\/l - k?), k) = sn(n, k) , (A.6) 



and has simple poles in u = 2nK(K) + (2m + 1) iK(\/l — k^), n,m S Z, with residue 
(— The functions cn and dn satisfy similar relations. Since am(ii, 0) = u, one has 
sn(n, 0) = sinn, cn(ti,0) = cosu and dn(ti,0) = 1. As k grows from to 1, the elliptic 
functions become more and more squarish. This is also apparent from their Fourier series, 
given by 

2KU) f2KU) \ 4^ q(2p+i)/2 
sn .- 1 — 



vr 



TT \ TT 



^ p=0 ^ 

2K(k) /2K(k) , \ 4^q(2p+i)/2 

cnf ^^V, = - E 1 + q2p+i ^°^((2y + ' 



2K(^ /2K(k) \ , , . 
dn( 1 = 1 + 42^ 2p Cos(2pV) , 



TT 



One should beware of the fact that some sources use m = k as parameter. 
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where q = q(K) is the elliptic nome defined by 



q = exp< — vr 



K(^ 



K(k) 



The elliptic nome has the asymptotic behaviour 



q(K) 



16 32 ^ ^ 



exp 



TT 



log[(l - k2)/16] 



l + O 



log^[(l-K2)/16] 



for K \ , 

for K y 1 . 



We also use the following identities, derived in |Jac691 p. 175]. For /c ^ 1, 



vr 



vr 



(A.8) 



(A.9) 



= C2A:,0 + ^C2fc,pY— ^COs(2pV) , (A. 10) 

^ p=l *^ 



where the C2kfl are positive constants (independent of ijj), and the other Fourier coefficients 
are given for the first few k by 



C2. 



C4, 



C6. 



(2p) , 



4 



vr 



3-4(2p)(l+^2^^ 2K(k) 
(2p)5- 20(2^)3(1 + ^2) 



(A.ll) 



2K(k)\^ , r, /2K(k'^'^' 

^ ^ ' +8(2p)(8 + 7k2 + 8k^ 



TT 



vr 



B Proofs of the Fixed-Point Argument 

In this appendix, we give the somewhat technical proofs of the fixed-point argument given 
in Section 13.61 We start by proving Lemma 13.61 stating a fixed-point equation equivalent 
to the stationarity conditions dGjy/dxq = 0, q £ Q. 

Proof of Lemma 13.61 The definitions (!3.54|) of an and /?„ imply, for any a,b^O, 



J2 J]p^^c^'?»(-+V2)ei'^og./M 



giv,9aeS«=i 
b 



E n 



Ltan(7rg^/iV) 



(B.l) 



It is more convenient to compute SGn / dp-^q rather than dGjy/dxq- We thus have to 
compute the derivatives of Qp with respect to p-g for all p. For p = 0, we have 



N 



£ 2^Go(^ + ^ "n) 



dp- 



n=l 



dP-q ' 



(B.2) 
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where (|B3]) shows that dun/dp-g = a;-9("+i/2) g-iV^og/M^ expand G'o(A + e'^On) into 
powers of e^, and plug in (|B.ip again. In the resulting expression, the sum over n vanishes 
unless (^i — (7 is a multiple of A^, say kN. This yields 

^ = e'^^^oA^/M ^ l^Gr^)(A)rK,(p) . (B.3) 

We consider the terms a = and a = 1 separately: 

• Since ^^^N+qiP) ~ ^kN-q vanishes for all k, the sum actually starts at a = 1. 

• The fact that T^^''^\p) vanishes whenever \£\ > N/2 implies that only the term k = 
contributes, and yields a contribution proportional to —pg. 

Shifting the summation index a by one unit, we get 



dp 



A similar computation for p 7^ shows that 



(B.4) 



Solving the stationarity condition 

= ^ = y e^iP^o ^ (B.6) 

with respect to pg, and singling out the term a = 6 = in (IB.SP to give the leading term 
then yields the result. □ 

The following estimates yield sufficient conditions for the operator T to be a contraction 
inside a certain ball, for the norm ||-||;^ introduced in (j3.72p . 

Proposition B.l. There exist numerical constants cq,ci > 0, such that for any A < Aq, 

and any N such that N e^^"^/'^^^ ^ 1/2, the estimates 



l + ^(llp||A+eAM?7(Ao,A))e||p|| 



W^P - VIU ^ ^S^^ \{\\ph V IIp'IIa) + sAMrjiXo, \)]e\\p - p'h (B.8) 



A|G^'(A)| 
\G'^{A)\ A4 

hold with ?7(Ao,A) = (e'^'/Ao) V (l/(Ao — A)), provided p and p' satisfy 

aV. ^0 - a a 



(B.7) 



^(MaVIIp'IU)^Co- lA^A-). (B.9) 



Proof: The lower bound 



ItanWAT)!^^^ 
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directly implies 



\ I 2M|p| 



e-^"|Pl . (B.ll) 



The assumption on N allows \pq\ to be bounded by a geometric series of ratio smaller than 
1/2, which is dominated by the term k = 0, yielding 

II (0) II / C2L0 (Ao-A)lgl/2M < C2L0 

where C2 > is a numerical constant. The fact that T^^'^\p) contains less than /v^'^+^'-i 
terms, together with (IB.lOh . implies the bound 

\rf^'\p)\ ^ AT'^+^'-ij^^y 6-^^/2^11^11^+" . (B.13) 

Assuming that \\p\\x ^ cpA^/Me for sufficiently small cq, it is straightforward to obtain 
the estimate 

In the sequel, we assume that ^ > 0, since by symmetry of the norm under permutation of 
Pq and p-q the same estimates will hold for q < 0. The norm of ^^"^^p, e) is more delicate 
to estimate. We start by writing 



where 

^«(^) = \\ + e-("°-')l^l «^p|-^(2M|p| + IfciV + g - 2Mp|)| . 

(B.16) 

We decompose 5'q(6) = S^{b) + 5"^ (6), where 5*+ (6) and 5^ (&) contain, respectively, the 
sum over positive and negative p. In the sequel, we shall only treat the term 5+ (6). The 
sum over k in (1B.16|1 is dominated by the term for which kN is the closest possible to 
2Mp — q, and can be bounded by a geometric series. The result for p > is 

Yexpi-^{2Mp + \kN + q - 2Mp\)\ ^ 04(6"^^ A e'^^/^M^ _ 

We now distinguish between two cases. 
• If q ^ 2M, we bound the sum over k by e"^^, yielding 

Stib) ^ ^— 2:/+^ e-^o^ ^ 4MC5^ . (B.18) 
Since e^^/^*^ ^ e^, it follows that 

|$(2)(. ^)\ < ^^^^^£2|| II -Ag/2Af 
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• If gr > 2M, we split the sum over p at q/2M. For 2Mp ^ q, we bound (1 + 2Mp/q) 
by 2 and the sum over k by q-M/'^m ^ p^j. 2Mp > g, we bound the the sum over k by 

St{b) ^ 2crM^^^^±i_e-^^/2M , (b.20) 

and thus 

Now (|R2T]) and ([Bl9]) . together with (jBrT2|) . imply ([B77|) . The proof of ([R8]) is similar, 
showing first the estimate 

-AE-=ik«l/2M| 



i=l i=l 



by induction on a, and then 



^ a(||/>|U V Wp'hr' ^-l..^.mi^M \\p _ p'll, (B.22) 



\Tt^'\p)-Tf''\p')\ ^ (a + fe)[iV(l|p||AV||p'|U)]'^+''-^(^y e-^'^l/2^||p-p'|U • (B.23) 

□ 



It is now easy to complete the proof of Proposition [ 
Proof of Proposition 13.81 Estimate ()B.7P for ||T p\\x implies that if 
, ^ ^0 . A3 / A|Gg(A)| \ 

then T{Bx{0,Rq)) C fiA(0,iio)- If in addition 

e^co— — lA " ^ A— , B.25 

then Estimate ()B.8|1 for ||Tp — T p'\\\ applies for p, p' G ;Ba(0, i?o)- It is then immediate to 
check that T is a contracting in B\{Q^Ro), a,s a consequence of ()B.24p . Thus the existence 
of a unique fixed point in that ball follows by Banach's contraction lemma. Finally, the 
assertions on the properties of p* follow from the facts that they are true for p^^\ that 
they are preserved by T and that p* = lim„^oo T^p^^^ ■ □ 
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